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Abstract 

This research deals with the mathematical modeling of the physical capital diffusion 
through the borders of the countries. The physical capital is considered an important 
variable for the economic growth of a country. Here we use an extension of the eco¬ 
nomic Solow model to describe how the smuggling affects the economic growth of the 
countries. In this study we rely on a production function that is non-concave instead 
of the classical Cobb-Douglas production function. In order to model the physical 
capital diffusion through the borders of the country, we developed a model based on a 
parabolic partial differential equation that describes the dynamics of physical capital 
and boundary conditions of Neumann type. Smuggling is present in many borders 
between countries and may include fuel, machinery and food. This smuggling through 
the borders is a problematic issue for the country’s economies. The smuggling problem 
usually is related mainly to a non-official exchange rate that is different than the official 
rate or subsides. Numerical simulations are obtained using an explicit finite difference 
scheme that shows how the physical capital diffusion through the border of the coun¬ 
tries. The study of physical capital is a paramount issue for the economic growth of 
many countries for the next years. The results show that the dynamics of the physical 
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capital when boundary conditions of Neumann type are different than zero differ from 
the classical economic behavior observed in the classical spatial Solow model without 
physical capital flux through the borders of countries. Finally, it can be concluded 
that avoiding the smuggling through the frontiers is an important factor that affects 
the economic growth of the countries. 

Keywords: Solow Model; Mathematical modeling; Physical Capital; Nu¬ 
merical Simulation; Neumann boundary conditions. 
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1 Introduction 

The Solow model presented in 1956 is the starting point for almost all analyses of economic 
growth. This dynamical model is still used in today’s macroeconomic theory despite its 
relative simplicity. Solow’s purpose was to develop a model to describe the dynamics of 
the growth process and the long-term evolution of the economy im H3- The Solow model 
describes the dynamics of the gross domestic product based on the variables labor, capital 
and technology. Thus, the classical Solow model can estimate the evolution of the physical 
capital over time through a differential equation. On the other hand, Solow model has been 
used in other areas such religion, ecommerce, internal brain drain, evaluation of the social 
benefit of water-saving agriculture and many others pi useei mg. 

The size and the rate of growth of the economy is of considerable interest for govern¬ 
ments and the public in general. The growth of the economy is related to many factors. 
One important factor for each country is the commercial relations with its neighbors. The 
smuggling phenomenon is a universal issue, which has been turned into a great problem in 
many countries. It is inappropriate since it has negative impacts on the economy, weakens 
the economical strategies and sometimes shows itself as a silent battle against the economy 
[6]. Porous borders and weak government capacity are important issues that are faced by 
several countries [2]. Furthermore, a large and growing underground economy may have 
serious economic implications on the countries economic growth. For instance, it can erode 
the tax base and cause important measurement errors in official economic statistics na- 
Because of its hidden nature, the underground economy is very difficult to measure. 

Yet from a much broader historical perspective, illicit cross-border flows of various sorts 
have been a defining feature of U.S. commercial relations with its neighbors from the very 
start, suggesting that there is much more continuity with the past than conventional accounts 
recognize. Porous borders and weak government capacity have long defined the region, and 
attempts to secure borders and tighten controls have often had the perverse and unintended 
consequence of creating a more formidable smuggling challenge. At the same time, efforts 
to regulate illicit border crossings have expanded the reach of central government authority 
and stimulated the development of border enforcement infrastructure and capacity. Bringing 
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this history back in to contemporary policy debates can offer fresh perspectives and lessons 
and provide an antidote to the often shrill and hyperbolic public discourse today about out 
of control borders. [21J 

On the other hand, petroleum product subsidies have increased in recent years. Many 
countries did not fully pass through the sharp increases in international petroleum product 
prices that occurred in 2007 and early 2008, resulting in a marked increase in subsidies. 
After declining along with oil prices during the second half of 2008, subsidies have again 
started to rise, renewing concerns about the fiscal costs. These concerns have been reinforced 
by the need in many countries to formulate an exit strategy from the recent crisis-related 
accumulation of public debt m- A subsidy by definition is any measure that keeps prices 
consumers pay for good or product below market levels for consumers or for producers above 
market. Subsidy means benefit given by the government to individuals or businesses whether 
in form of cash, tax reduction or by reducing the cost of goods and services. The purpose of 
subsidy is to help individuals and businesses purchase/acquire essential goods and services 
that they may not be able to afford, under normal circumstances. Subsidies take different 
forms and some have a direct impact on price. These include grants, tax reductions and 
exemptions or price controls ESI- Fuel subsidy is applied in oil producing countries, such 
as: Venezuela, Iran, Saudi Arabia, Egypt, Burma, Malaysia, Kuwait, China, Taiwan, South 
Korea, Trinidad and Tobago, and Brunei |2j2j. 

Generally, there are two main government subsidies: fuel subsidy where consumers pay a 
fraction of the price that consumers are supposed to pay and price controls regarding several 
important products including food. Subsidy, in economic sense, exists when consumers of 
a given commodity are assisted by the government to pay less than the prevailing market 
price of same. [23J 

Based on the previous aspects some countries such as Venezuela are facing an extraor¬ 
dinary smuggling of subsidized petroleum products and other regulated priced products 
through the borders. This economic behavior is usual in countries with lower prices in 
comparison with their border countries. The countries with higher prices are often substan¬ 
tial beneficiaries through cross-border smuggling. When relevant, governments should also 
highlight that subsidies promote smuggling, shortages and black markets dzi As has been 
mentioned, the supply and demand for smuggled goods depend on interregional price dispar¬ 
ities in the presence of a trade ban H. The smuggling through the borders of Venezuela has 
been increased in the last few years. People from Venezuela, Colombia and Brazil often par¬ 
ticipate in the informal circulation of legal and illegal goods across the border. This economic 
behavior is also presented in other places. Due to the discrepancy between local understand¬ 
ing of legality and national laws, many potentially newsworthy illegal exchanges are not 
addressed in the media [20]. There are transborder roads between Colombia and Venezuela. 
However, the main one connects the Venezuelan city of San Antonio with Cucuta, the capi¬ 
tal of the Colombian province Norte de Santander. The Simon Bolivar international bridge 
can be crossed by foot, car or bicycle. Moreover, no documents are required and there are 
few inspections through the bridge. On the Colombian side, small-scale vendors of gasoline, 
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known as pimpineros , peddle the fuel along the highway, without safety precautions, from 
small plastic containers. A gallon (3.8 litres) of gasoline costs the equivalent of 4.92 dollars 
in Colombia, compared to just eight cents in Venezuela, because of the heavy subsidies that 
make gasoline in this country among the cheapest in the world [29]. The price differential 
allows 1,000 percent profits to be made from smuggling gasoline. Moreover, nowadays it is 
not only gasoline, but also food and medicines, subsidized in Venezuela, that are trafficked. 

Fuel smuggling to Colombia costs Venezuela US 1.4 billion per year, with 30,000 barrels 
smuggled daily, according to a report coming out of a bilateral meeting between the Colom¬ 
bian and Venezuelan governments in August, 2013 [3TJ. In addition, fuel smuggling to other 
countries such Brazil, and the islands of Curazao, Aruba, and Trinidad and Tobago are also 
made. Recently, international studies points to the existence of a correlation between glob¬ 
alization and trans-border economic crimes PI- Moreover, some authors have argued that 
the human use of the single international political economy, which globalization signifies its 
transition, has receded to the logic of thinking globally and acting locally [ 16j . 

The purpose of this paper is to shed some new light on the debate regarding smuggling 
through the borders of the countries. In this paper we study the effects of the smuggling 
of goods across the borders of countries. We rely on an extension of the economic Solow 
model to describe how the smuggling affects the economic growth of the countries. The 
model is based on a parabolic partial differential equation that describes the dynamics of 
physical capital and the boundary conditions are of Neumann type in order to model the 
physical capital diffusion through the borders of the country. However, there exist many 
other options that has been proposed to deal with this type of economic issues and is still 
an open economic debate ra- 

2 Spatial Solow model 

The Solow model describes the dynamics of the gross domestic product based on the vari¬ 
ables labor, capital and technology. Thus, the classical Solow model can estimate the time 
evolution of the physical capital through a differential equation. In regard to physical capital 
there is not complete agreement about what this term refers and how to measure it [!24j. In 
fact some authors have introduced the concept of social capital 1233 However, in general 
for economist, physical capital refers to a factor of production (or input into the process of 
production), such as machinery, buildings, or computers [28] . In economic theory, physical 
capital is one of the three primary factors of production, also known as inputs production 
function. Another capitals considered in economy are human capital which is the result of 
investment in the human agent and financial capital. Human capital is the stock of compe¬ 
tencies, knowledge, social and personality attributes, including creativity, embodied in the 
ability to perform labor so as to produce economic value [28]. 

As has been stated in many articles, the mathematical expression for the total production 
or output in terms of the input variables labor L(t), physical capital K{t) and technology 
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A(t) is given by, 


Y(t) = A(t) F(K(t), 


(1) 


where F denotes the production function. On the other hand, the stock of capital depreciates 
over time at a constant rate S and only a fraction of the output ( cY(t ) with 0 < c < 1) is 
consumed and sY(t) is the savings (s=l-c). Thus, the Solow model describes the change of 
the physical capital as, 

'FFl = S Y(t) - mt). ( 2 ) 

In addition, using the production function one gets, 


'FFl = sA{() F{K{()< £(()) _ sK(t). 


(3) 


Now if we introduce a change of variable k(t) = j4J in order to compute the per capita 
physical capital one obtains, 


mt) = F(K(t)/L(t), 1). (4) 

Thus, one gets that, 

< ^^- = sA(t)f(k(t))-(S + n(t))k(t), (5) 

where n(t) denotes the labor growth rate. Notice that the labour-augmenting technology or 
technology progress A(t) affects positively the variation of physical capital and, depreciation 
and labor growth rate affects negatively. 

There are several variations that may be included for the based Solow model ft For 
instance if we consider a net flow of physical capital to a given location or space interval 
we need to extend the Solow model in order to account for the flow ME]. Thus, k(x, t) 
denotes the physical capital stock held by a representative household located at x at the 
time t, in a bounded domain f2 G R n . Moreover, here we consider that the production 
function f(k(t )) is the same whatever is the location x. In addition, the initial physical 
capital distribution, k(x, 0), is assumed to be known and C°. Thus, one gets the following 
model for the dynamics of the physical capital, 


- ^ ^ — A k(x, t ) = sA(x, t ) f(k(t )) — 5k(x, t), 0 x [0, T] 

(J L 

S7k(x,t ) = h(k(x,t)), dfl x [0,T] 

k(x, 0) = ko(x), x G fl. (6) 

This model has been called a spatial Solow model 0E]. Notice that here we introduce a 
Neumann type boundary conditions in order to deal with the flow of physical capital through 
the borders of the country in an illegally way. Moreover, the boundary conditions are related 
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in dQ with the initial physical capital distribution in order to have a well posedness model. In 
addition, the model (|6]) is still open since the production function f(k(t)), the technological 
progress function A(x,t ) and the initial physical capital distribution k 0 (x) are unknown. 

Several economic growth models are based on production functions. These functions play 
a crucial role in the economy and different types of function have been used. It specifies the 
maximum output for all possible combinations of input factors and therefore determines the 
way the economic model evolves in time [B EDI IS]. The Cobb-Douglas production function is 
especially notable for being the first time an aggregate or economy-wide production function 
had been developed, estimated, and then presented to the profession for analysis; it marked 
a landmark change in how economists approached macroeconomics ra Furthermore, the 
function was not developed on the basis of any knowledge of engineering, technology, or 
management of the production process. Nowadays, the Cobb-Douglas production function 
is by far the most used to represent the technological relationship between the amounts of 
two or more inputs, particularly physical capital and labor, and the amount of output that 
can be produced by those inputs. The Cobb-Douglas production function is non-negative, 
increasing and concave, and, satisfies the well-known Inada conditions mm- However, 
other productions functions may be used in conjunction with economic models such the 
Solow model. For instance Leontief production or nonconcave functions are accepted in the 
applied literature. Here we rely on a particular production function introduced in |9j. This 
function is called an S-shaped production function (nonconcave) allowing the Solow model 
to have a richer dynamics which includes the existence of a poverty trap. As it has been 
mentioned in [9], the complexity of the Solow model using the nonconcave function suggest 
to do the analysis by means of numerical techniques. Thus, for our numerical simulations 
we consider the following nonconcave production function 

ak p 

1 TW q ' 

where the parameters a, (3, p and q are usually determined by adjusting the production 
function by least squares to economic data. In regard to the technological progress function 
A(x,t ) many mathematical expressions have been used. However, a constant value has been 
assumed in several works. For instance in [8] authors used a function such the economy has 
a technological center, that may coincide with the geographical center. On the other hand, 
in [H] it has been assumed that the technological progress can be modeled by a diffusion 
process coupled to the physical capital diffusion. Thus, these authors assume that A(x, t ) is 
an unknown to be determined. In our numerical simulations of the next section we assume 
for the function A(x,t ) different mathematical expressions including constant in space and 
time, and varying linearly or exponentially with time as has been suggested in the literature. 

There are lots of evidence regarding smuggling activities in many countries however, there 
are no reliable statistics available which can be quite misleading |6|. Numerical simulations 
are performed in next section in order to study different economic scenarios. As it has 
been mentioned the parameter values are difficult to estimate for each country and general 
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numerical values are assumed without loss of generality. In regard to the initial physical 
capital distribution k 0 (x), we use in the numerical simulations four different types. The first 
one is a constant but is only consistent when the function h(x,t ) involved in the Neumann 
type boundary conditions is equal to zero. The other three types are consistent with h(k, x ) 7^ 
0. The first one has a Gaussian form with higher values of physical capital in the middle 
of the physical space. The last two functions are given by parts, where one is composed by 
linear forms and the last one is more complex since the function h(k, x) depends explicitly 
on the physical capital available on the border. More details regarding these mathematical 
expressions are given in the next section. 

3 Numerical simulations 

We perform several numerical simulations of the spatial Solow model for different scenarios 
in order to study the effect of physical capital diffusion through the borders of the countries. 
The proposed model is based on a parabolic partial differential equation and the boundary 
conditions are of Neumann type. The numerical results are obtained using an explicit finite 
difference scheme and small time steps in order to fulfill the Courant Friedrich Lewy condition 
for the numerical stability of the solution. However, implicit methods without stability 
problems to numerically solve the spatial Solow model can be found in [zuHuaj. 

We study four important factors of the spatial Solow economic model such the depreci¬ 
ation rate, technological progress, initial conditions and boundary conditions. For all these 
factors we use different numerical values or type of functions in order to observe the effect 
of each one on the physical capital dynamics. We simulate several unit times in order to 
observe transient and steady states. It is important to remark that in this first study we 
do not rely on real world data since our first main goal here is to develop the methodology 
and analyze the effect of physical capital diffusion through the borders of the countries. 
Moreover, real world numerical values are very difficult to estimate and for some variables 
such the level of technology is not observable and estimated values need to be used [4]. On 
the other hand, estimation for the production function also would be necessary and a huge 
amount of economic panel data is required which is out of the scope of this study. 

Here we show the behavior of the physical capital under different scenarios regarding the 
boundary conditions which are related to the borders of the countries. For each scenario 
numerical simulations are performed in order to understand the effect of each boundary 
condition on the dynamics of the physical capital. For our numerical results we consider 
the S-shaped nonconcave production function 0 with a = 0.0005, (3 = 0.0005, p — 4 and 
q = 4. Finally, it is important to mention that a spatial dimension of L = 100 is used for 
all the numerical simulations without loss of generality and a numerical value for the saving 
rate s = 1 is taking for simplicity on the analysis of the dynamics. It is well-known that the 
saving rate s in the long term has no effect on the growth rate of the economy [9]. 


3.1 No physical capital flux through the borders 

In this first case we consider the classical spatial Solow model without flux on the borders of 
the countries. In addition, we consider here that the initial condition for the physical capital 
is given by a uniform distributed function k(x, 0) = 100 and that no technological progress 
is present, i.e. A(x,t) = 1. The boundary conditions in explicit form for the physical capital 
model (|6j) are given by, 


dk(0 , t) 


= 0 


and 


dk(L,t) 


= 0 . 


( 8 ) 


dx dx 

In these first numerical simulations we would like to point out the effect that the depreci¬ 
ation rate S has on the physical capital dynamics. It can be seen in Fig. [I] that the physical 
capital decreases to zero for depreciation rates of 6 = 0.05 and 6 = 0.5. However, for the 
latter case the rate of decrease is much higher due to the depreciation. It is important to 
remark that one way to avoid the physical capital decrease is to increase the technological 
progress A(x,t ) of the country. 



Time t 



Time t Space x 


Figure 1: Numerical solutions of the spatial Solow model without physical capital flux 
through the borders of the country. Initial condition for the physical capita is given by 
a uniform distributed function k(x, 0) = 100, no technological progress is included here, i.e. 
A(x,t ) = 1. The increase of the depreciation rate 8 implies a faster decrease of the physical 
capital as it can be seen on the left (5 = 0.05) and right (5 = 0.5) hand-sides. 


3.2 No physical capital flux through the borders with technolog¬ 
ical progress 

In this second case we again consider the classical spatial Solow model without flux on the 
borders of the countries. However, in this case we consider that the technological progress 
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is increasing every year (time). In regard to initial condition we set an uniform distributed 
function k(x, 0) = 100. For the technological progress we consider two scenarios; the first 
one with a linear increase A(x,t ) = t and the second one with an exponential increase of the 
technological progress A(x , t) = e°' 01t . In these scenarios we would like to point out the effect 
that the technological progress A(x,t ) has on the physical capital dynamics. It can be seen in 
Fig. [2] that in both scenarios the physical capital decreases initially and then after some time 
the physical capital starts to increase due to the technological progress. Notice that in the 
previous cases, where no technological progress was included, the physical capital drops to 
zero after some time. Thus, it is important to remark that the technological progress A(x,t ) 
of the country can help for the country’s growth, since this variable is related to physical 
capital [H EDI 0 m El- Furthermore, it has been shown that spatial diffusion mechanisms 
of technology progress allows the develop of other regions and its depends on the spatial 
distance JHUH2]. Technological progress may be obtained by investment in human resource, 
specialization and industrial structure. 



Figure 2 : Numerical solutions of the spatial Solow model without physical capital flux 
through the borders of the country. Initial condition for the physical capita is given by 
a uniform distributed function k(x, 0 ) = 100 , no technological progress is included here, 
i.e. A(x,t ) = e 0 01 b The increase of the depreciation rate 6 implies a faster decrease of the 
physical capital as it can be seen the left (S = 0.05) and right (S = 0.5) hand-sides. 


3.3 Physical capital flux through the borders. 

In this case we would like to study the effect that has the physical capital flux through the 
borders using the spatial Solow model. In this scenario we consider that the technological 
progress is constant in order to focus in the effect that has the capital flux through the 
borders. We focus in three types of boundary conditions in order to simulate the flow 
of physical capital and such they are consistent with h(k(x,t )) 7 ^ 0. The first one has a 
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Gaussian form with higher values of physical capital in the middle of the physical space in 
order to simulate a country where usually the regions near the borders are more exposed to 
the consequences of the smuggling. For the first case we consider the following boundary 
conditions, 


dk(0 , t) 
dx 


D 0 


and 


dk(L,t) 

dx 


(9) 


where the parameter D 0 = lOe - ^ 502 )/ 1000 measures the amount or intensity of physical capital 
flow through the borders of the country. In regard to initial condition we set a Gaussian 
distributed function k(x, 0) = 100 e - G- 50 ) / d ; where d is a parameter that measures the 
differences between the central region and the border regions, where a large value of the 
parameter d means small differences for the region’s physical capital. Thus, initial physical 
capital decreases as we depart from the center. It is important to notice that initial condition 
must be different than a constant uniform distribution like in the previous cases in order 
to have a consistency with the physical capital flux through the borders. In other words, 
the flux induces a nonuniform distribution of the physical capital. Here we consider also 
two scenarios with different deprecation rates 5 = 0.05 and 8 = 0.0005. It can be seen in 
Fig. 1 that when a depreciation rate is high the physical capital will decrease steadily to 
zero. However, for a lower value of the depreciation rate the physical capital will increase 
despite the flux of it through the borders. However, due to the physical capital flux now a 
lower depreciation rate is required in order to avoid a decrease for the capital. Thus, it is 
important to remark that the flux of physical capital through the borders of the country can 
undermine the country’s growth. 
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Figure 3: Numerical solutions of the spatial Solow model with physical capital flux through 
the borders of the country. Initial condition for the physical capital is given by a Gaussian 
distributed function k(x, 0) = 100 e - ^ -50 ) 2 / 1000 , A(x, t) = 1 and 8 = 0.05, 0.0005 respectively. 


In the second scenario with physical capital flux we consider that the initial condition is 
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given by the following function given by parts, 


k(x, 0) 


mx if 0 < x < 10 

C max if 10 < x < 90 

—m(x — 90) + C max if 90 < x < 100. 


where we assume C max = 1000 and m = 100. It is important to notice that this physical 
capital initial condition varies with the position on the regions near the borders and in the 
central regions is assumed constant in order to simulate as close as possible a real scenario. 
Moreover, this initial condition is consistent with the boundary conditions that allow the flux 
of physical capital through the borders. For this case we consider the following boundary 
conditions, 


dk( 0, t) 


= m 


and 


dk(L,t) 


—m, 


( 10 ) 


dx dx 

where the parameter m measures the amount or intensity of physical capital flow through 
the borders of the country. In order to study the effect of this type of initial condition 
we also include the effect of the deprecation rate with two numerical values S = 0.05 and 
5 = 0.005. It can be seen in Fig. [4] that for the largest depreciation rate the physical capital 
will decrease steadily to zero. On the other hand, for the smallest value of the depreciation 
rate the physical capital increases despite the flux of it through the borders. Notice that the 
physical capital is larger in the central regions than in the border regions due to the flux on 
the borders. 



Figure 4: Numerical solutions of the spatial Solow model with physical capital flux through 
the borders of the country. Initial condition for the physical capita is given by a function 
given by parts, A(x,t ) = 1 and S = 0.05,0.005 respectively. 

In our last scenario we introduce a more complex boundary condition in order to avoid 
possible negatives values for the physical capital at the long term and to analyze the ef¬ 
fects that has the physical capital flux, depreciation and the initial condition on the econ- 
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omy(physical capital). In this scenario we consider that the initial condition is the following 
function given by parts, 


k(x, 0) 


pDx 

e D( 100—x) 


if 0 < x < 10 
if 10 < x < 90 
if 90 < x < 100. 


where we set the parameter D = 1. This initial condition is consistent with the boundary 
conditions that allow the flux of physical capital through the borders. For this case we 
consider the following boundary conditions, 


dk( 0, t) 
dx 


Dk( 0, t) 


and 


dk(L,t) 

dx 


—Dk( 0, t), 


( 11 ) 


where the parameter D measures the fraction of physical capital that flow through the 
borders of the country. 

In order to study the effect of the deprecation rate with take two numerical values 8 = 0.05 
and 8 = 0.005. It can be seen in Fig. [5] that for the largest depreciation rate the physical 
capital decreases in with a lower rate in comparison with the case for the smallest value. 
Notice again that the physical capital is larger in the central regions than in the border 
regions due to the flux on the borders. In this way, as has been mentioned previously 
the flux of physical capital through the borders is an economical issue that can affects the 
economic growth of the countries and its growth. 



Figure 5: Numerical solutions of the spatial Solow model with physical capital flux through 
the borders of the country. Initial condition for the physical capita is given by a function 
given by parts, A(x,t ) = 1 and 8 = 0.05,0.005 respectively. 
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4 Conclusions 

In this study we proposed specific boundary conditions of Neumann type for the spatial 
Solow model to explain and understand the effect of physical capital diffusion through the 
borders of the countries. The physical capital has been considered an important variable 
for the economic growth of a country. In the last decade fuel, medicine, machinery and 
food smuggling through the borders of the country has been a problematic issue for the 
Venezuelan economy. Moreover, in the last 12 months the smuggling problem in Venezuela 
has been increasing dramatically due mainly to a non-official exchange rate that is now more 
than ten times the official rate. 

Here we have used an extension of the economic Solow model to describe how the smug¬ 
gling of fuel, machinery, goods and others elements affects the economic growth of the coun¬ 
tries. In this study we have relied on a production function that is non-concave instead of the 
classical Cobb-Douglas production function. The proposed model was based on a parabolic 
partial differential equation that describes the dynamics of physical capital and the boundary 
conditions were of Neumann type in order to model the physical capital diffusion through 
the borders of a country. The numerical results were obtained using an explicit finite dif¬ 
ference scheme that shows how the diffusion of physical capital to the border countries is a 
paramount issue for the economic growth of the countries. 

We can summarize the results using four important factors such the depreciation rate, 
technological progress, initial conditions and boundary conditions. We found that when 
there is not out flux of physical capital the country would need less effort regarding the 
variables depreciation rate, technological progress to achieve an increase of capital. On the 
other hand, when flux of capital is considered thorough the borders of the country, it is more 
demanding to increase the physical capital since a lower depreciation rate and an increasing 
function for the technological progress may be necessary for the economic growth of the 
country. 

Finally, we found that initial condition for the distribution of the physical capital in 
conjunction with the boundary condition have some effect on the dynamics of the physical 
capital on the long term. This fact differs with the scenario when there is not flux since at 
the long term the physical flux is distributed uniformly. Thus, it is important to remark that 
the flux of physical capital through the borders of the country can undermine the country’s 
growth. Furthermore, numerical simulations showed that the physical capital is larger in 
the central regions than in the border regions due to the flux on the borders. In this way, 
we can conclude that the flux of physical capital through the borders is an economical issue 
that may affect Venezuela’s economy and its growth. 

Further research in this economic issue should lead to improve the understanding of the 
economic problem that many countries face. Moreover our results can help government eco¬ 
nomic institutions and population to deal and propose economic policies to tackle the flux of 
physical capital through the borders of the country such as smuggling of fuel, machinery and 
food. The model also gives likely future outcomes as well as the implications of alternative 
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economic policies. A natural continuation of this research is the extension to the Ramsey 
model and the introduction of different production functions. 


A Appendix 


In order to compute the numerical simulations of problem @ , we define the domain dis¬ 
cretization as Qq = [0, L] x [0,T], We discretize the domain O 0 in the following form: we 
chose positive integers N, T 0 such that the spatial step size is given by Ax = L/N, while the 
time step size is At = T/T 0l and we put Xj = jAx, t n = nAt, 0 < j < N, 0 < n < T 0 . 
We seek approximations of the solution at these mesh points, these approximate values will 
be denoted by k™ ~ k(xj,t n ). By replacing the first derivative and the second order space 
derivative by the classical difference quotient it follows that 


dk(xj,t n ) _ k(xj,t n+ 1 ) - k(xj,t n 


+ O(At), 


dt At 

d^kjx^tn) _ k{xj +1 ,t n ) - 2 k(xj, tn) + k{Xj-i,t n ) 

dx 2 Ax 2 

For the sake of simplicity, we define the difference operator 

h n — Oh n h n 
K j +1 AK 3 + K i -1 


+ 0(Ax 2 




Ax 2 


1 <j<N. 


( 12 ) 


(13) 


Using (12) and (13), it obtains an explicit standard numerical scheme for the problem 


Q, such that for 0 < n < T 0 , 0 < j < N , 

k™ +1 = (1 - A tS)k] + AtSj + AtsA] + f(k]), Q x [0, T], 
kj = ko(xj), Xj e U. 


(14) 


Now, to calculate the values of the operator 5” in (13) for j = 0, and j = N requires 
us to invoke the fictitious values and k r ^ +l . These approximations are determined using 
the boundaries conditions, such that for this case are: 


dk{ 0, t) 
dx 


= dohi(k(0,t)), 


dk(L, t) 
dx 


= d L h 2 (k(L,t)). 


(15) 


kx +1 , the following expressions 


Next, approximating the equations (15) with central differences scheme, one gets for /e” 1 and 

(16) 


k-i = K + 2d 0 Axh 1 (ko), k^ +1 = k^-i - 2d L Axh 2 (k^). 
Thus, we calculate in the boundary points using 

k™ +1 = (1 - A t8)kfi + 2At — - + ^ xhl ^ k - + A tsA$ + f{k$), 


(17) 
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_i_i k^i, i — bn + di,Axh2{k%) 

k n N +1 = (1 - At5)k n N + 2A+ /(*&), 0 < n < T 0 . 

Thus, we obtain an explicit numerical scheme which is well-known to be numerically stable 
and convergent whenever The numerical errors are proportional to the time step 

and the square of the space step. 
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